International Journal of Modern Physics B, 
© World Scientific Publishing Company 



BETHE ANSATZ AND 
CLASSICAL HIROTA EQUATION 



P. WIEGMANN 
James Franck Institute and Enrico Fermi Institute 
University of Chicago, Chicago, IL 60637, USA 

We discuss an interrelation between quantum integrable models and classical soliton 
equations with discretized time. It appeared that spectral characteristics of quantum in- 
tegrable systems may be obtained from entirely classical set up. Namely, the eigenvalues 
of the quantum transfer matrix and the scattering S-matrix itself are identified with a 
certain r-functions of the discrete Liouville equation. The Bethe ansatz equations are ob- 
tained as dynamics of zeros. For comparison we also present the Bethe ansatz equations 
for elliptic solutions of the classical discrete Sine-Gordon equation. The paper is baned 
on the recent study of classical integrable structures in quantum integrable systems. El 

1. Introduction 

1. In 1981 Hirotai proposed a difference equation which unifies aU known continuous 
soHton equations. A particular case of the Hirota equation is a bihnear difference 
equation for a function r(n, Z,rn) of three discrete variables: 

ar(n, I + 1, m)T(n, /, m + 1) + (3T{n, /, m)T{n, I + l,m + 1) 

+ 'jT{n + 1,1 + 1, m)T{n - 1 J, m + 1) = , (1) 

where it is assumed that a+/3+7 — 0. Different continuum limits at different bound- 
ary conditions then reproduce continuous soliton equations (Kadomtsev-Petvi^shvili 
equation, Toda lattice, etc.). On the other hand, T{n,l,m) can be identified^l with 
the r-function of a continuous hierarchy. 

The same equation (with a particular boundary condition) has suddenly ap- 
peared in the theory of quantum integrable systems as a -fusion relation for the 
transfer matrix (trace of the quantum monodromy matrix) .ffl 

Thus classical integrability emerges in quantum systems not as a limiting proce- 
dure but rather as inherent and exact. It appears that while solving an integrable 
quantum problem we in fact deal with classical integrable equations, but with a 
discretized time. Indeed, the Bethe ansatz equations and the eigenvalues of the 
S-matrix and transfer matrix may be obtained by solving classical Hirota equation. 

At present, we cannot treat an appearance of the classical (discrete time) equa- 
tions in quantum problem better than an observation, although all this indicates 
that classical and quantum integrable systems have much more in common than a 
regular classical limit Ti ^ 0. 

2. The transfer matrix is one of the key objects in the theory of quantum inte- 
grable systemsJj Transfer matrices form a commutative family of operators acting 
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in the Hilbert space of a quantum problem. Let Ri^^iu) be the i?-matrix acting in 
the tensor product of Hilbert spaces Vi (8) Va- Then the transfer matrix is a trace 
over the auxiliary space of the monodromy matrix that is the matrix product 
of N i?-matrices with the common auxiliary space: 

TA{u\yi) = RN,A{u-yN)---R2,A{u-y2)Ri,A{u-yi), 

Ta{u) = i^ATA{u\yi). (2) 

The transfer matrices commute for all values of the spectral parameter u and dif- 
ferent auxiliary spaces: 

[Ta{u),Ta'{u')]^0. (3) 

They can be diagonalized simultaneously, i.e. their eigenstates do not depend on 
A and u. The family of eigenvalues of the transfer matrix is the first question to 
be addressed for an integrable system - the spectrum of a quantum problem can be 
expressed in terms of eigenvalues of the transfer matrix. 

Moreover, the transfer matrix may be used to obtain the scattering S-matrix. 
Indeed, at u = the fundamental i?-matrix becomes a permutation R(u = 0) = P. 
Therefore, choosing A to be the same as the representation of , say, the first particle 
and setting u = yi, the transfer matrix becomes the S'-matrix of scattering the first 
particle with rapidity j/i by the other — 1 particles 

SA{u\y2, --yN) = Rn,a{u - yN) ■ ■ ■ R2,a{u ~ y2) ■ (4) 

3. The transfer matrix corresponding to a given representation in the auxiliary 
space can be constructed out of transfer matrices for some elementary space by 
means of the fusion procedurexi'QU The fusion procedure is based on the fact that 
at certain values of the spectral parameter u the i?-matrix becomes essentially a 
projector onto an irreducible representation space. The fusion rules are especially 
simple in the ^i-case. Consider for example the rational solution of the Yang-Baxter 
equation for the fundamental (spinor) representations of SU(2): 

Ris{u) = u-2P (5) 

where P is the permutation operator. At a special value of the spectral parameter 
u = ±2, the i?-matrix becomes a projector onto irreducible moduli of the tensor 
product [1/2] ® [1/2] = [0] © [1] - the singlet (spin-0) at and onto the triplet (spin-1) 
:i?(±2) — 4PT. Then the transfer matrix ^^^^(u) with spin-1 auxiliary space is ob- 
tained from the product of two spin-1/2 monodromy matrices 3f (u) with arguments 
shifted by 2: 

T^{u) = tr[i](i?i,i(-2)fi(i/ + l)fi(7. - l)i?i,i(-2)) . 

A combination of the fusion procedure and the Yang-Baxter equation results in 
numerous functional relations (fusion rules) for the transfer matrixE They were re- 
cently combined into a universal bilinear form.clEl The bilinear functional relations 
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have the most simple closed form for the models of the A^^-i-series and representa- 
tions corresponding to rectangular Young diagrams. 

Let r° (u) be the transfer matrix for the rectangular Young diagram of length a 
and height s. They obey the following bilinear functional relation: 

T^{u + l)T^{u - 1) - T^+,{u)T!_,{u) = T^+\u)Tr\u) . (6) 

Since Tg{u) commute at different u, a, s, the same equation holds for eigenvalues 
of the transfer matrices, so we can (and will) treat T^{u) in as number- valued 
functions. The bilinear fusion relations for models related to other Dynkin graphs 
were suggested in ref.tl. 

Remarkably, the bilinear fusion relations appear to be identical to the Hirota 
equation (]l|). Indeed, one can eliminate the constants a, /?, 7 by the transformation 

(-a/7)"'/2 

so that 

T„(/ + l,m)Tn{l,m + 1) - T„(/, m)Tn{l + l,m + 1) 

- T„+i(Z + l,m)T„_i(/,m + 1) = (7) 

and then change variables from light-cone coordinates n, m to the direct variables 

a = n, s = I + m, u = I — m ~ n, 

Tn{l,m)=Ti\^{l-Tn-n). (8) 

At least at a formal level, this transformation provides the equivalence between (^) , 
(0) and §). 

Relation between Hirota's equations and fusion relations goes much further. It 
is known that a general class of Hirota's equations are related with the geometr 
Grassmannian manifolds and can be constructed out of a general Young tableau.1 
If, for example, the Young diagrams Y = [a^^ , (01-1-02)''^] consist of two rectangular 
blocks (one with ai lines of length si + S2 and the second with 02 lines of length 
Si), then the higher Hirota equations hold 

= Kri:s'{u~i)T:;L7X\^+^)+K\:^h7 (9) 

This is also a fusion relation for the transfer matrix of the auxiliary space of the 
representation Y. 

Leaving aside more fundamental aspects of this "coincidence" , we exploit, as a 
first step, some technical advantages it offers. Specifically, we treat the functional 
relation (|^) not as an identity but as a fundamental equation which (together with 



4 P. Wiegmann 



particular boundary and analytical conditions) completely determines all the eigen- 
values of the transfer matrix. The solution of the Hirota equation then appears in 
the form of the Bethe ansatz equations. We anticipate that this approach makes it 
possible to use some specific tools of classical integrability and, in particular, the 
finite gap integration technique. 

The origin of Tg{u) as an eigenvalue of the transfer matrix imposes specific 
boundary conditions: 

T^^O, a/ 0,1,.., A: (10) 

and, what is equally important, certain analytical properties. As a general con- 
sequence of the Yang-Baxter equation, the transfer matrices may be always nor- 
malized to be elliptic polynomials in the spectral parameter, i.e. finite products of 
Weierstrass cr-functions (as in (^9|) below). Solutions with that kind of analytical 
properties (called elliptic sokdunpaLare well-known in the theory of classical soliton 
equations due to the works. Ej't^llSlla 

4. To illustrate elliptic solutions in classical continuum time soliton equations 
let us recall the result of the Refills for the Kadomtsev-Petviashvili (KP) equation. 
A particular elliptic solution can be written in the form 

N 

u{x,y,t)^2Y,V{x^x,{y,t)) (11) 

i=l 

where V is Weierstrass elliptic function and all dynamics in t and y is hidden in the 
behaviour of poles Xi{y,t). A very existence of this ansatz is an implementation of 
the integrability. It appears that poles obey the elliptic Calogero-Moser many-body 
system with respect to y 

x': = -j2p'ix.,-x,). (12) 

i 

Elliptic solutions can be lifted up to the discrete caseEB Then poles Xi (or roots 
of the elliptic polynomial) obey a fully discretized version of the Calogero-Moser 
model, which can be recognized as the Bethe ansatz equations. 

In this paper we shall consider the discrete version of the classical Liouville 
equation. Its elliptic solutions gives the eigenvalues of the transfer matrix and 
the 5-matrix of the quantum integrable systems associated with Ai algebra. For 
extension to Ak-i algebra see Ref.tl. We also compare it, but do nrit study in detail, 
with Bethe ansatz solution of the discrete Sine Gordon cquation.tZl 

The Hirota equation for the Ai case may be obtained from the more general (^), 
by setting the boundary conditions 

r; = 0, a ^0,1, 2 (13) 

Then T^'^ obey the discrete Laplace equation while the T^{u) = Ts{u) obeys the 
discrete Liouville equation (we continue to refer it as the Hirota equation): 

T,{u + l)T,{u - 1) - T,+i{u)T,^i{u) = T°iu)T^{u). (14) 
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In fact the values of (u) and T^{u) can be easily determined. They are the 
transfer matrices in most symmetric and antisymmetric states of the A2 case. These 
states have dimension 1. In the most symmetrical state there is no scattering, so 
Tg{u) = 1. The value of the transfer matrix in the most antisymmetric state is 
known as the quantum determinantEa 

-i-A^) — —/ r—/ 7T- (15) 



where 

7(w- 



7(u) 

and 



(16) 



N qk-1 

</-(«) =n n <^iv{u-yk+i)). (17) 

k=l i=-(jfc + l 

Here cr(x) is the Weierstrass a-function and is a spin of the fcth quantum space. 
The roots yk and degree TV of the elliptic polynomial </<(«) and parameters of the 
Weierstrass function are parameters of the quantum model. Below we concentrate 
mostly on the simplest case gfc = 1. 

Apparently, different solutions of have a number common zeros and (u) at 
general parameters yk are polynomials of degree s'^qk- These zeros can be gauged 
out, so all Ts{u) (different eigenvalues at different s and qk) remain polynomials of 
the same degree N. The maximal gauge transformation is 

s N 9fc-s 

T4u) ^ T4u) Y[^{u + s-2p)Y[ Yl a{T]{u - s + qk + 2 - 2p) 

p=2 k=l p=l 

where the last product holds for q^ > s. This gauge transformation simplifies the 
r.h.s. of (|l|) 

Ts{u + \)Ts{u - 1) - Ts+i{u)Ts^i{u) = (t>{u + s)(f){u ~ s - 2). (18) 

In this normalization all are polynomials of the same degree with no common 
zeros. 

A general solution of the discrete Liouville equation ( p^ is parametrized by two 
arbitrary functions of one variable. Not all of them correspond to eigenvalues of 
the quantum transfer matrix, but only elliptic solution. Elliptic polynomials may 

is) 

be characterized by its roots 

N 

where Ag, /i^ do not depend on u. Similar to the continuum case, roots obey some 
dynamics in a discrete "time" variable s. It is determined by the Bethe- ansatz-like 
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equations 

^(4'^ + + 1) ^ fr a(ry(zf - 4-^^ - 1)) 

^{^? -s + l) _ A,^, ^ 'yMzf ^ + 1)) 

0(z} ^ - s - 1) ^s+l cr(?7(2] ^ - ^ - 1)) 

The constant Ag can be easily found in the rational and trigonometric case by 
comparing the leading powers of polynomials in the (^8|). In the rational case it 
obeys the equation 

Al - As+iA,^i = 0. 

Each solution of these equations gives an eigenvalue of the transfer matrix. Equa- 
tions ( |20| ) and (|2l| ) determine the zeros of the transfer matrix. The traditional 
Bethe ansatz equations 



e 



(Uj) -r-r (T{r]{Uj - Itfc + 2)) 



4>{uj - 2) o-(»7("i - Uk - 2)) 

corresponding to zeros of the Baxter Q matrix are discussed in Sec. 5. 
2. Other Forms of the Hirota Equation 

The equation ( p^ ) is known as a discrete version of the Liouvillc equation^ written 
in terms of the r-function. It can be recast to a more universal form in terms of 
the discrete Liouvillc field 

' ^ ' (j){u + s)<f){u - s - 2) ^ ' 



which hides the function <f){u) in the r.h.s. of ([181). The equation becomes 

Y,iu - l)Y,iu + 1) = (r,+i(M) + l)(i;-i(w) + 1) . (24) 

with a boundary condition Yo{u) — 0. (Let us note that the same functional equa- 
tion but with different analytic properties of the solutions appears in the thermo- 
dynamic Bethe ansatz.0) 

The functional equations ( p^ ) can be further rewritten in an integral form. In- 
troducing Ss(u) — \iiYs{u), one gets 

/ioG 
9^rn.^(.,^.A + + es-i)) = E Ss.^^'^iu) (25) 

where (pk{u) are some known functions computed out of the r.h.s. of the (p^, such 

as (pkiu + 1) + (fkiu -1) = 0n 
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3. Difference and Continuum ^i-Equation 

Liouville equation. In the continuum limit one should put Ys{u) — S~'^e~'^^^'*\ 
u = 5~^x, s = S^^t. An expansion in 5 — > then gives the continuous Liouville 
equation 

92^-92^ = 2exp((p). (26) 

Sine- Gordon equation. For illustrative purposes let us compare it with discrete 
Sine-Gordon (SG) equation, which requires quasi-periodic boundary condition with 
respect to the Dynkin diag rami 

Tf+i(u) = e"A2°Tf-i(u - 2), (27) 

where a and A are parameters. Plugging this condition into (^, we get: 

T^iu + l)T}iu - 1) - Tl^,{u)Tl,iu) = e-X'T^iu)T^iu - 2), (28) 

T^{u + l)T^{u - 1) - T^^,{u)Tl,iu) e--T}iu)TXu + 2). (29) 
Let us introduce two fields p^'" and on the square (s, u) lattice 

T^iu) = exp(p^'" + ^^'-), (30) 

T}{u + l) = X'/^eMp''^-^'n, (31) 

and substitute them into ( p8|) and (p9|). Finally, eliminating p*'", one gets the 
discrete SG equation: 

sinh ((^''+i'" + <y5'*-i^" - (/j-^'^+i - (^^'"-1) 

= A sinh ((^''+1'" + + (p^^'^+i + if'-''-^ + a) . (32) 

The constant a can be removed by the redefinition — > yj*'" — a/4. 
Another useful form of the discrete SG equation appears in variables 



Under condition (^^ one has 

X!+\u) = Xri(u - 2), A2x;+1(u + 1)X,"(7.) = 1 , (34) 
so there is only one independent function 

Xl{u) = x,{u) = -e-°'X-^ exp ( - 2(p"'" - 2ip'''''-^) . (35) 
The discrete SG equation becomesB^Elii 

, . , . (A + x,(u+l))(A + x,(u-l)) 
x.+r (u)x._, (u) = + • (36) 
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In the limit A ^ ( |3q ) turns into the discrete LiouviUe equation ( p4| ) for Ys{u) = 
— 1 — X~^Xs{u). 

4. Determinant representations 

Relation between the Hirota equation and Plucker relations of the coordinates of 
the Grassmannian manifolds suggests^ a numerous determinant representations of 
its solutions. The most familiar onecil allows one to express T^{u) through T^{u) 
or T^{u). For instance the determinant formula giving the evolution in a holds 



T^iu) = det {T^+'-'{u + 



(37) 



Another representation appeared to be more suitable for particularities of the 
LiouviUe equation (|^) and (p^. Xtey explicitly express its solution in terms 
of 2k arbitrary functions of one variablea as a determinant of the k x k matrix 



Ta{l,m) 



detM,% 



hi{u + s + a + 2j) if j = 1, k — a; i = 1, k 
hi{u — ,s + a + 2j) if J = fc — a + 1, fc; i ~ 1, . 



In the Ai case the determinant has a particular simple form 
T°(w) = <i>{u + s) 



R{u + s) Q{u + s) 
R{u + 2 + s) Q{u + 2 + s) 



T^{u) = ^iu-s-2) 



R{u - s - 1) Q{u - s - 1) 
R{u - s + l) Q{u -s + 1) 



Ts{u) 



Q{u + s + l) R{u + s + l) 
Q{u — s) R{u — s) 



Compare first two equations one finds 

g(u) = g(w~l), R{u)^R{u~l) 

so that 



Ts{u) 



Q{u + S + 1) R{u + s + 1) 
g(w-s-l) R{u-s~l) 



(38) 

(39) 

(40) 
(41) 

(42) 
(43) 



5. Bethe-Ansatz 

The functions Q,R are further determined by the r.h.s. of the ( p^ ) and by the 
analyticity requirement that Tgiu) is an elliptic polynomial (p^). One can show 
that it also means that Q{u) or R{u) (but not necessarily all together) is an elliptic 
polynomials of degree N/2 multiplied by exponential function 

M 

Q{u)=e'^^^\{a{i^{u-u,)) (44) 
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Let Uj, j — 1,...,M be zeros of Q{u). Then, evaluating ( |39| ) aX u — uj, 
u = Uj — 2, we obtain the relations 

= Q{uj + 2)R{uj) , (j){uj - 2) = -Q{uj - 2)R{u.j) , (45) 

whence it holds 

(/.(w, -2) Q(w, -2)' ^ ^ 

Equation ( ^6| ) is the celebrated Bethe ansatz equation (p2|). 

In the elliptic case the degrees of the elliptic polynomial Q{u) (for even N) is 
equal to M = N /2^ provided that rj is incommensurable with the lattice spanned by 
wi, UJ2- In the trigonometric and rational cases there are no such strong restrictions 
on degrees M and M of Q- This is because a part of their zeros may tend to 
infinity thus reducing the degree. Whence M and M can be arbitrary integers not 
exceeding TV. However, they must be complement to each other: M ~\- M — N . 
The traditional choice is M < N/2. In particular, the solution Q{u) = 1 {M = 0) 
corresponds to the simplest reference state ( "bare vacuum" ) of the model. 

6. Linear Problem 

The function Q naturally appears as the object of the linear problem of the discrete 



Liouville equation. Indeed, the equation (14) admits a zero curvature representa- 



tion. Let us start from a more general equation 

Ts{u + l)Ts{u - 1) - Ts+i{u)Ts-i{u) = (j){u + s)0(u - s), (47) 

with two independent functions ((>, (j) and later impose the relation (/'(u) = (l){u — 2). 
The nonlinear equation appears as a consistency conditions of the following auxiliary 
linear problems 

Ts+i{u)Q{u + s) - Ts{u - l)Q{u + s + 2) = 0(m + s)Q{u - s - 1), (48) 

T,+i{u)Q{u - s + 1) - Ts{u + l)Q{u - s - 1) = 0(u - s)Q{u + s + 2) , (49) 

where we introduced two functions of one variable Q{u) and Q{u). These are the 
functions appeared in the determinant representations (^9[), (|40|), and (^l]). 

The fact that functions of the linear problem depend on one variable is a specific 
feature of the Liouville equation. 

To compare let us present the linear problem for the discrete Sine-Gordon equa- 
tion (H) and (PI). It reads 



TI+,{u)F!{u) - Tl{u - l)FO+i(t. + 1) = T^An)Fl+,{u), (50) 

TiVi (^ - 2)Fl {u) - T^iu ~ 3)i^,Vi + 1) = Tl {u)F^^, {u - 2) (51) 

T^^,{u)F^{u) - T^u - 1)F!+,{u - 1) = e-'^T}{u)F}^,iu + 2), (52) 

T,Vi(^ - l)F}{u) - T}iu)F}+,{u - 1) = e'^X^T^iu - 3)F^+,iu) (53) 
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In this case the wave functions F^(u) and F}{u) depend on both s and u. 
7. The Baxter Relation 

A general solution of the discrete Liouvillc equation (for arbitrary (j) and (f>) may be 
expressed through two independent functions Q{u) and Q{u). Below we follow the 
sam e li n e s developed for solving continuous classical Liouville equation (see, e.g., 
Refs£3'E3 and references therein). First we rearrange (^) and ( p9| ) as 

(t){u - 2)Q{u + 2) + (l>{u)Q{u - 2) = A{u)Q{u), (54) 

4>iu)Q{u + 3) + 4>{u + 2)Q{u - 1) = A{u)Q{u + 1), (55) 

where we have introduced the quantities 

(j){u - 2)Ts+i{u - s) + (j){u)Ts-i{u ~ s - 2) 

= — , (56) 

ls(u - s - 1) 

^ ^ 4>{u + 2)Ts+i{u + s) + ^{u)Ts-i{u + S + 2) 



Due to the consistency condition (47), A{u) and A{u) are functions of one variable 
and do not depend on s. The symmetry between u and s allows one to construct 
similar objects which in their turn do not depend on u. Functions A{u) and A[u)^ 
in the r.h.s. of ( p^ ) and ( |55| ) are the conservation laws of the s-dynamics. 

Running ahead, let us note that the connection between and cf), 4>{u) — <j){u—2), 
and its consequence T_i(ii) = 0, simplifies ( |5^ ) - (|5^). Putting s — and using the 
boundary condition T_i(u) = 0, we find 

A{u) = A{u) = Ti{u) . (58) 

Therefore, we have the relations 

Ts{u - l)Ti{u + s) = 0(u + s - 2)T,+i(u) + ,^(u + s)r,_i(u - 2), (59) 

Ts{u + l)ri(M - s) = (t){u - s)Ts+i{u) + (t>{u-s^ 2)Ts-i{u + 2), (60) 

0(u - 2)Q{u + 2) + <j!)(u)Q(ii - 2) = Ti{u)Q{u) . (61) 

These en uaJities are also known as fusion relations and have been obtained in 
RefsE'EZlta, and ( |6l| ) is also known as the Baxter T-Q relation.t^ So the Baxter 
Q function and T-Q relation naturally appear in the context of the auxiliary linear 
problems for the bilinear fusion relation. 

The Baxter equation ( |6l| ) gives rise to a generalized spectral problem for a 
discrete operator of the second order: find all elliptic "eigenvalues" of the operator 
(f){u — 2)e^^" + 0(u))e~^^". One may show that in this case the eigenfunction Q{u) 
itself is an elliptic polynomial. The Bethe ansatz is a practical tool to solve this 
problem. 
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Let US consider ( |54| ) (resp. (p5D) as a second order linear difference equation, 
wliere ttie function A{u) {A{u)) is determined from the initial data. Let R{u) (resp. 
R{u)) be a second (linearly independent) solution of ( p^ (resp. (p5|)) normalized 
so that the wronskians are 



W{u) = 
W{u) = 



R{u) 

R{u) 
R{u + '. 



Q{u) 
Q{u + \ 

Q{u) 
Q{u + 2) 



= (t>{u), 
<j){u + ; 



(62) 
(63) 



Then the general solution of (^) is given by the determinant (^Tj). 

For any given Q{u) and Q{u) the second solution R{u) and R{u) (defined modulo 
a linear transformation R{u) R{u) + aQ{u) ) can be explicitly found out of the 
first order recurrence relations (p9|), (pO[), if necessary. Then (^Tj) determines Ts{u) 
only in terms of Q{u): 



(64) 



This formula has been obtained in RefsEJ'Ea by direct resolving fusion recurrence 
relations (|5^) and (pO|). The Bethe ansatz equations may be also obtained as a 
conditions of cancelation of poles in (|6^. 

Let us list some more useful representations. 

The conserved quantities A{u) and A{u) in terms of Q and R have the form: 



A{u) = Q{u + 2)R{u - 2) - R{u + 2)Q{u - 2), 

A{u) = R{u + 3)Q{u - 1) - Q{u + 3)R{u - 1), 

which are direct corollaries of (^4|) and (|55|). 

In its turn Baxter's relation (|6l|) and its "chiral" versions (| 
enjoy some determinant representation^ 



Ts{u) 




Ts+i{u - 


1) 


Q{u + 


s + 


1) 




Ts+iiu + 


1) 






Q{u + 


s + 


3) 




Ts+2 ( u + 


2) 


Ts+3{u + 


1) 


Q{u + 


s + 


5) 




Ts{u) 




Ts+i{u + 


1) 


Q{u~ 


s) 






Ts+i{u - 


1) 


Ts+2{u) 




Q{u- 


s — 


2) 




Ts+2 [u - 


2) 


Ts+siu - 


1) 


Q{u^ 


s — 


4) 





(65) 
(66) 

and (|55|) also 



(67) 



(68) 



provided that Ts{u) obeys the Hirota equation. 
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8. Zeros of the transfer matrix 

The linear problem (^) and ( ^ ) are convenient to find the Bethe-ansatz-like equa- 
tions 

® for zeros of Ts{u). Substituting u = z^'^'^^^ and u = -f- 1 into ( p8| ) and 
u — zf*"''^'' and u = z^^'''' — 1 into ( ^ ) we obtained the ( pO| ) and ( [2l| ) referred in the 
introduction. In the rational case (cr(?7w) —^u)As = s + l. 

Let us compare the Bethe-ansatz-like equations of the Liouville equation ( pO| ) 
and ( pl| ) with the Sine-Gordon equation ( p8| ) and (|29|). Let ^^'(s) and zj{s) be the 
zeros of the elliptic polynomial 

N 

Substituting u = ^^'^(s) into the second pair of equations of the linear problem 
and ( p3| ) and excluding F's after some algebra one obtains 



n 



L V^(rykHs + 1) - z°{s + 2) - 3); V^(^kHs + 1) - z}{s + 1) + 1) 
a(r7(zKs + 1) - + 2) - 2) _ A^^A]^, 



N 

n 



a(r7(zO(s + 1) - z^s + 1) + 2)\ /airj{z^is) - zO(,s + 1) + 1) 



(70) 



a(r;(zO(s + 1) - zi(s + 1) + 1)7 V f^('7(^."(5) ~ ^"(s) + 1) 
ct(77(zO(,s + 1) - zjjs + 2) - 1) ^ AgAi+i 
^ a(,y(zO(. + 2)-zK^ + 2)) ^2+1^^ 



(71) 



The constants As may be easily found in the rational or trigonometric limits by 
comparing leading orders of polynomials in the Hirota equation. In the rational 
case one gets: 

{Air-Al^,AU^e-X'{A'ir (72) 

{A'sr-A's+,AU=e-'^>^'{Air. (73) 

These equations gives a discrete dynamics of zeros in s. They are to be 
with dynamics of zeros of elliptic solutions of classical nonlinear equations 



9. Crossing Symmetry and the S-Matrix 

As it was noticed in the introduction, once the eigenvalues of the transfer matrix are 
found it is easy to get eigenvalues of the S'-matrix (Q). In relativistic field theories 
the S'-matrix is normalized in a way to be unitary and crossinvariant, i.e. to be 
unitary in both s and t channels. 

Let us consider for instance the spin 1/2 fundamental Ai rational S'-matrix. It 
has the form: 

S{u) = U{u)R{u) (74) 



Bethe Ansatz and Classical Hirota Equation 13 

where the i?-matrix is given by (^) and the ampHtude U(u) is restricted by t he 
unitarity and the cross symmetry. The unitarity conditions in both channels ara23 

Uiu)U*i~u) - ! — (75) 

(p(u + 2)(f)(-u + 2) 

Uiu + 2) ^ ^{u + 2)<^(u - 2) 

U{u-2) (/.(m)0(u + 4) ■ ^ ' 

The minimal solution is 

n(u) = '^^ 4 ''-y 4 ) 



The Hirota equation has a natural normalization which gives the cross- 
invariant unitary S'-matrices in all representations. It is a gauge which brings the 
quantum determinant (|l^) to 1 and the Hirota equation and the Baxter relation 
( |6l|) to the universal form 

Ts{u + 1)ts{u - 1) - Ts+i{u)ts-i{u) = 1 (78) 

Q{u + 2) + Q{u-2)=Ti{u)Q{u). (79) 
This is achieved by the gauge 

T.(u)=r.(.):g^ (80) 

where 

fiu+l)fiu-l)^jiu)j{u + 4). (81) 

Then Ts is a meromorphic function with prescribed positions of poles. It brings the 
i?-matrix to the cross unitary ^-matrix 



This factor can be identified with the factor U in (|25|). Indeed, (81) implies 

f±fl^^(.)0(. + 4) 
f{u + 3)f{u-3) ^iu)<f>iu + A) 



fiu+l)fiu-l) ^{u-2)cf>iu + 2) 
This is equivalent to the unitarity and cross- unitarity relation (^5|) and ([76|). 

10. Conclusion and Outlook 

It turned out that quantum integrable models have a deep connection with classical 
discrete soliton equations: the fusion rules for quantum transfer matrices are iden- 
tical to the Hirota bilinear difference equation with a certain boundary conditions 



14 P. Wiegmann 



and analytical requirement. This coincidence goes far beyond the simplest example 
considered in this paper. It most likely holds for all Lie algebras and representations 
and all integrable boundary conditions. 

Under this identification eigenvalues of the transfer matrix are represented as 
T-functions. Positions of zeroes of the solution is determined by the Bethe ansatz 
equations. We therefore conclude that the information one usually obtains from 
the Yang-Baxter equation may be also found from classical discrete time equations, 
without implying a quantization. 

In general the Bethe ansatz technique, what has been thought of as a specific tool 
of quantum integrability is shown to exist in classical discrete nonlinear integrable 
equations. We demonstrated this on two examples - discrete Liouville and discrete 
Sine-Gordon equation. It is tempting to compare solutions of the quantum Sine- 
Gordon field theory and discrete classical Sine-Gordon equation. In a forthcoming 
paper we extend the theory of elliptic solution known in continuum soliton equations 
to a more general class of discrete equations by means of the Bethe ansatz. 

The difference bilinear equation, although with different analytical requirements, 
has appeared in quantum integrable systems in another context. Spin-spin corre- 
lation functions of the Ising model obey^hilinear difference equation that can be 
recast into the form of Hirota equation.ElSEj More recently, nonlinear equations for 
correlation functions have been derived for rnore general class of quantum integrable 
models, by virtue a new approach of the RefJ23._Tliermodynamic_Bethe ansatz equa- 
tions written in a form of functional relationsc3 (see, e.g., Ref.c3) appeared to be 
identical with Hirota equation (although with another analytic properties) as well. 

All these suggest that the Hirota equation may play a role of a master equa- 
tion simultaneously for both classical and quantum integrable systems, such as the 
"equivalence" between quantum systems and discrete classical dynamics may be 
extended beyond the spectral properties discussed in this paper. 
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